Sexually-dimorphic targeting of functionally-related genes in COPD

Background: There is growing evidence that many diseases develop, progress, and respond to therapy differently in men and women. This variability may manifest as a result of sex-specific structures in gene regulatory networks that influence how those networks operate. However, there are few methods to identify and characterize differences in network structure, slowing progress in understanding mechanisms driving sexual dimorphism. Results: Here we apply an integrative network inference method, PANDA (Passing Attributes between Networks for Data Assimilation), to model sex-specific networks in blood and sputum samples from subjects with Chronic Obstructive Pulmonary Disease (COPD). We used a jack-knifing approach to build an ensemble of likely networks for each sex. By adapting statistical methods to compare these network ensembles, we were able to identify strong differential-targeting patterns associated with functionally-related sets of genes, including those involved in mitochondrial function and energy metabolism. Network analysis also identified several potential sex- and disease-specific transcriptional regulators of these pathways. Conclusions: Network analysis yielded insight into potential mechanisms driving sexual dimorphism in COPD that were not evident from gene expression analysis alone. We believe our ensemble approach to network analysis provides a principled way to capture sex-specific regulatory relationships and could be applied to identify differences in gene regulatory patterns in a wide variety of diseases and contexts. Electronic supplementary material The online version of this article (doi:10.1186/s12918-014-0118-y) contains supplementary material, which is available to authorized users

Epidemiology, genetics, and subtyping of preserved ratio impaired spirometry (PRISm) in COPDGene.

BackgroundPreserved Ratio Impaired Spirometry (PRISm), defined as a reduced FEV1 in the setting of a preserved FEV1/FVC ratio, is highly prevalent and is associated with increased respiratory symptoms, systemic inflammation, and mortality. Studies investigating quantitative chest tomographic features, genetic associations, and subtypes in PRISm subjects have not been reported.MethodsData from current and former smokers enrolled in COPDGene (n = 10,192), an observational, cross-sectional study which recruited subjects aged 45-80 with ≥10 pack years of smoking, were analyzed. To identify epidemiological and radiographic predictors of PRISm, we performed univariate and multivariate analyses comparing PRISm subjects both to control subjects with normal spirometry and to subjects with COPD. To investigate common genetic predictors of PRISm, we performed a genome-wide association study (GWAS). To explore potential subgroups within PRISm, we performed unsupervised k-means clustering.ResultsThe prevalence of PRISm in COPDGene is 12.3%. Increased dyspnea, reduced 6-minute walk distance, increased percent emphysema and decreased total lung capacity, as well as increased segmental bronchial wall area percentage were significant predictors (p-value <0.05) of PRISm status when compared to control subjects in multivariate models. Although no common genetic variants were identified on GWAS testing, a significant association with Klinefelter's syndrome (47XXY) was observed (p-value < 0.001). Subgroups identified through k-means clustering include a putative "COPD-subtype", "Restrictive-subtype", and a highly symptomatic "Metabolic-subtype".ConclusionsPRISm subjects are clinically and genetically heterogeneous. Future investigations into the pathophysiological mechanisms behind and potential treatment options for subgroups within PRISm are warranted.Trial registrationClinicaltrials.gov Identifier: NCT000608764

‘Location, Location, Location’: a spatial approach for rare variant analysis and an application to a study on non-syndromic cleft lip with or without cleft palate

Motivation: For the analysis of rare variants in sequence data, numerous approaches have been suggested. Fixed and flexible threshold approaches collapse the rare variant information of a genomic region into a test statistic with reduced dimensionality. Alternatively, the rare variant information can be combined in statistical frameworks that are based on suitable regression models, machine learning, etc. Although the existing approaches provide powerful tests that can incorporate information on allele frequencies and prior biological knowledge, differences in the spatial clustering of rare variants between cases and controls cannot be incorporated. Based on the assumption that deleterious variants and protective variants cluster or occur in different parts of the genomic region of interest, we propose a testing strategy for rare variants that builds on spatial cluster methodology and that guides the identification of the biological relevant segments of the region. Our approach does not require any assumption about the directions of the genetic effects. Results: In simulation studies, we assess the power of the clustering approach and compare it with existing methodology. Our simulation results suggest that the clustering approach for rare variants is well powered, even in situations that are ideal for standard methods. The efficiency of our spatial clustering approach is not affected by the presence of rare variants that have opposite effect size directions. An application to a sequencing study for non-syndromic cleft lip with or without cleft palate (NSCL/P) demonstrates its practical relevance. The proposed testing strategy is applied to a genomic region on chromosome 15q13.3 that was implicated in NSCL/P etiology in a previous genome-wide association study, and its results are compared with standard approaches. Availability: Source code and documentation for the implementation in R will be provided online. Currently, the R-implementation only supports genotype data. We currently are working on an extension for VCF files. Contact: [email protected]

WISARD: workbench for integrated superfast association studies for related datasets

Background: A Mendelian transmission produces phenotypic and genetic relatedness between family members, giving family-based analytical methods an important role in genetic epidemiological studies—from heritability estimations to genetic association analyses. With the advance in genotyping technologies, whole-genome sequence data can be utilized for genetic epidemiological studies, and family-based samples may become more useful for detecting de novo mutations. However, genetic analyses employing family-based samples usually suffer from the complexity of the computational/statistical algorithms, and certain types of family designs, such as incorporating data from extended families, have rarely been used. Results: We present a Workbench for Integrated Superfast Association studies for Related Data (WISARD) programmed in C/C++. WISARD enables the fast and a comprehensive analysis of SNP-chip and next-generation sequencing data on extended families, with applications from designing genetic studies to summarizing analysis results. In addition, WISARD can automatically be run in a fully multithreaded manner, and the integration of R software for visualization makes it more accessible to non-experts. Conclusions: Comparison with existing toolsets showed that WISARD is computationally suitable for integrated analysis of related subjects, and demonstrated that WISARD outperforms existing toolsets. WISARD has also been successfully utilized to analyze the large-scale massive sequencing dataset of chronic obstructive pulmonary disease data (COPD), and we identified multiple genes associated with COPD, which demonstrates its practical value. Electronic supplementary material The online version of this article (10.1186/s12920-018-0345-y) contains supplementary material, which is available to authorized users

Heme metabolism genes Downregulated in COPD Cachexia.

IntroductionCachexia contributes to increased mortality and reduced quality of life in Chronic Obstructive Pulmonary Disease (COPD) and may be associated with underlying gene expression changes. Our goal was to identify differential gene expression signatures associated with COPD cachexia in current and former smokers.MethodsWe analyzed whole-blood gene expression data from participants with COPD in a discovery cohort (COPDGene, N = 400) and assessed replication (ECLIPSE, N = 114). To approximate the consensus definition using available criteria, cachexia was defined as weight-loss > 5% in the past 12 months or low body mass index (BMI) (< 20 kg/m2) and 1/3 criteria: decreased muscle strength (six-minute walk distance < 350 m), anemia (hemoglobin < 12 g/dl), and low fat-free mass index (FFMI) (< 15 kg/m2 among women and < 17 kg/m2 among men) in COPDGene. In ECLIPSE, cachexia was defined as weight-loss > 5% in the past 12 months or low BMI and 3/5 criteria: decreased muscle strength, anorexia, abnormal biochemistry (anemia or high c-reactive protein (> 5 mg/l)), fatigue, and low FFMI. Differential gene expression was assessed between cachectic and non-cachectic subjects, adjusting for age, sex, white blood cell counts, and technical covariates. Gene set enrichment analysis was performed using MSigDB.ResultsThe prevalence of COPD cachexia was 13.7% in COPDGene and 7.9% in ECLIPSE. Fourteen genes were differentially downregulated in cachectic versus non-cachectic COPD patients in COPDGene (FDR < 0.05) and ECLIPSE (FDR < 0.05).DiscussionSeveral replicated genes regulating heme metabolism were downregulated among participants with COPD cachexia. Impaired heme biosynthesis may contribute to cachexia development through free-iron buildup and oxidative tissue damage

COPD association and repeatability of blood biomarkers in the ECLIPSE cohort.

BACKGROUND: There is a need for biomarkers to better characterise individuals with COPD and to aid with the development of therapeutic interventions. A panel of putative blood biomarkers was assessed in a subgroup of the Evaluation of COPD Longitudinally to Identify Surrogate Endpoints (ECLIPSE) cohort. METHODS: Thirty-four blood biomarkers were assessed in 201 subjects with COPD, 37 ex-smoker controls with normal lung function and 37 healthy non-smokers selected from the ECLIPSE cohort. Biomarker repeatability was assessed using baseline and 3-month samples. Intergroup comparisons were made using analysis of variance, repeatability was assessed through Bland-Altman plots, and correlations between biomarkers and clinical characteristics were assessed using Spearman correlation coefficients. RESULTS: Fifteen biomarkers were significantly different in individuals with COPD when compared to former or non-smoker controls. Some biomarkers, including tumor necrosis factor-α and interferon-γ, were measurable in only a minority of subjects whilst others such as C-reactive protein showed wide variability over the 3-month replication period. Fibrinogen was the most repeatable biomarker and exhibited a weak correlation with 6-minute walk distance, exacerbation rate, BODE index and MRC dyspnoea score in COPD subjects. 33% (66/201) of the COPD subjects reported at least 1 exacerbation over the 3 month study with 18% (36/201) reporting the exacerbation within 30 days of the 3-month visit. CRP, fibrinogen interleukin-6 and surfactant protein-D were significantly elevated in those COPD subjects with exacerbations within 30 days of the 3-month visit compared with those individuals that did not exacerbate or whose exacerbations had resolved. CONCLUSIONS: Only a few of the biomarkers assessed may be useful in diagnosis or management of COPD where the diagnosis is based on airflow obstruction (GOLD). Further analysis of more promising biomarkers may reveal utility in subsets of patients. Fibrinogen in particular has emerged as a potentially useful biomarker from this cohort and requires further investigation. TRIAL REGISTRATION: SCO104960, clinicaltrials.gov identifier NCT00292552.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Analyzing networks of phenotypes in complex diseases: methodology and applications in COPD

Background: The investigation of complex disease heterogeneity has been challenging. Here, we introduce a network-based approach, using partial correlations, that analyzes the relationships among multiple disease-related phenotypes. Results: We applied this method to two large, well-characterized studies of chronic obstructive pulmonary disease (COPD). We also examined the associations between these COPD phenotypic networks and other factors, including case-control status, disease severity, and genetic variants. Using these phenotypic networks, we have detected novel relationships between phenotypes that would not have been observed using traditional epidemiological approaches. Conclusion: Phenotypic network analysis of complex diseases could provide novel insights into disease susceptibility, disease severity, and genetic mechanisms